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The hypothesis that dark matter and dark energy are unified through the Chaplygin gas is 
reexamined. Using a generalization of the spherical model which incorporates effects of the 
acoustic horizon we show that an initially perturbative Chaplygin gas evolves into a mixed 
system containing cold dark matter in the form of gravitational condensate. Furthermore, by 
including both condensate and residual gas, we demonstrate that the observed CMB angular 
and baryonic power spectra are reproduced 



An appealing scenario in which dark matter and dark energy are different manifestations of 
a common structure , may be realized through the Chaplygin gas, an exotic fluid obeying 

P = -A/p , (1) 

which has been extensively studied for its mathematical properties^. The cosmological potential 
of Eq. was first noted by Kamenshchik et aJ^ who observed that the Chaplygin gas interpo- 
lates between matter with p ~ a -3 at high redshift and a cosmological constan t like p ~ \f~A~ as 
a tends to infinity. Of particular interest is that Eq. (^Q) may be obtained from^^^ 



c m = -Va Ji - e,„ e, v , (2) 



by evaluating the stress-energy tensor T^ u , introducing = 9, a / \l g a @ 0, a 9,g for the four- 



velocity, and p = v A/ y/1 — g^ u 9,^ 9, u for the energy density. One recognizes £bi a Born-Infeld 



type Lagrangian familiar in the Z)-brane constructions of string/M theory^. The Lagrangian 
(j2J is a special case of the string-theory inspired tachyon Lagrangian^ in which the constant 
\/~A is replaced by a potential V(0). 

To be able to claim that the Chaplygin gas (or any other candidate) actually achieves unifi- 
cation, one must be assured that initial perturbations can evolve into a deeply nonlinear regime 
to form a gravitational condensate of superparticles that can play the role of cold dark ma tter 
(CDM). In comoving coordinates, the solution for inhomogeneous Chaplygin gas cosmology iS^ 



A + Bh. (3) 

Here 7 is the determinant of the induced metric jij = giogjo/goo — 9iji an d B can be taken as 
constant on the scales of interest. Eq. © allows us to implement the Zel'dovich approximation^ 
the transformation from Euler to Lagrange (comoving) coordinates induces jij = dkiD^Dj" , 
where ZV = a{5i — bip/) is the deformation tensor, <p is the velocity potential, and the quantity 
b = b(t) describes the evolution of the perturbation. The Zel'dovich approximation offers a means 
of extrapolation into the nonperturbative regime with the help of Eq. (j3J) and 

7 = a 6 (l - Ax6) 2 (l - A 2 6) 2 (l - A36) 2 , (4) 

where the Aj are the eigenvalues of (p^ 3 . When one (or more) of the A's is (are) positive, a 
caustic forms on which 7 — > and p/p — > 0, i.e., at the locations where structure forms the 
Chaplygin gas behaves as dark matter. Conversely, when all of the A's are negative, a void 
forms, p is driven to its limiting value VA, and the Chaplygin gas behaves as dark energy, 
driving accelerated expansion. For the issue at hand, the Zel'dovich approximation has the 
shortcoming that the effects of finite sound speed are neglected. Indeed, in the Newtonian limit 
p <C p, an explicit solution for the perturbative density contrast of the pure Chaplygin gas 

5 peit (k,a) oc a~ 1/4 J 5 / M (d s k) , (5) 

has been obtained^. Here J u (z) is the Bessel function, k the comoving wave number, and d s 
the comoving sonic horizon given by 

d - 2 ( 1 -^ 1/2 q7/2 (6) 



with the equivalent matter fraction = y/B/(A + B) = VB/p cr . Thus, for d s k <C 1, <5 pe rt ~ a, 
but for d s k S> 1, <5 pe rt undergoes damped oscillations. 

Since the structure formation occurs in the decelerating phase, we can address the question 
within Newtonian theory by generalizing the spherical model. In the case of vanishing shear 
and rotation, the continuity and Euler-Poisson equations become 

p + 3Hp = 0; 3W + 3W 2 + 4vrGp + V- ^Vp^J =0, (7) 

where TL is the local Hubble parameter. It is reasonable to approximate S(t, x) = p(t, x)/p(t) — 1 
by the spherical lump 5(t, x) = 6k(t) sin(kx)/(kx) . We then find that ^( a ) satisfies^ 

2x " , 3 , 3 X 4 (a5' k ) 2 49 / a \ 7 5 k 



where = (d s k) 2 1 1 a. Eq. © reproduces © at linear order and extends the spherical dust 
model by incorporating the Jeans length through the last term. 
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Figure 1: (a) Top left: Evolution of 6k (a) in the spherical model, Eq. ®, from a cq = 3 x 1CT 4 for a k = 0.0476, 
<5fc(ffleq) =0.004 (solid) and <5t,(a eq ) =0.005 (dashed), (b) Bottom left: The critical <5fc(a oq ) versus a k . (c) Top 
right: CMB angular power spectrum for a mixture of Chaplygin gas and condensate with f2 gas = 0.7 and h — 0.7, 
varying Ocdm and v 2 . (c) Bottom right: Baryon power spectra for the mixed Chaplygin gas with Ocdm = 0.2 
and v 2 — 0.9995 and for the ACDM model. A normalization factor of 1.35 is used. 



In Fig. ^ we show the evolution of two initial perturbations from radiation-matter equality 
for ak = a rc ion = 1/21. One sees that for sufficiently small 5k(a eq ), the acoustic horizon can 
stop 5k (a) from growing even in a mildly nonlinear regime. Conversely, for 5k(a cq ) above the 
threshold, 5k(a) — > oo at finite a just as in the dust model. The critical 5k(a eq ) dividing the two 
regimes depends strongly on (Fig. \J3p). Qualitatively similar conclusions have been reached 
in a different way by Avelino et Since the critical 5k i s co mmensurate with the peak in the 
conditional probability distribution for spheroidal collapse^^, and ak is only weakly dependent 
on the comoving wave number, we can thus be confident that the Chaplygin gas will evolve 
at high redshift into a mixed system consisting of smoothly distributed gas and gravitational 
condensate. The latter will participate in hierarchical clustering as CDM. 

Homogeneous world models containing a mixture of c old dark matter and Chaplygin gas 
have be en successfully confronted with lensing statisticsU^ as well as with supernova and other 
tests^. Clearly, it is only a matter of interpretation to replace 'cold dark matter' by 'Chaplygin 
droplet matter' (condensate). On the other hand, "unification" h as ofte n been misconstrued to 
mean a mixture of baryonic and perturbative Chaplygin gas only^^^l Owing to the damped 
oscillation of perturbations and the driving decay of the gravitational potential, it is hardly 
surprising that the baryon plus the purely perturbative Chaplygin gas model is in gross conflict 



with the CMB^^^and mass power spectrum!^ data. Hence we have undertaken a calculation 
of the CMB anisotropics and the mass power spectrum in our unification scenario based on the 
Chaplygin gas and including nonlinear condensate. As the sonic horizon is negligible at recombi- 
nation and the droplets affect the CMB and the power spectrum only through the gravitational 
potential, it is adequate to treat them throughout as ordinary cold dark matter, parametrized 
by f^cDMj in the perturbation equations. The residual uncondensed gas is parametrized as 



Pgas(a) = Pcr^gas V 2 + (1 - V 2 ) /a 6 ' (9) 



in the background and 



v z (a)=v 2 v z + (l-v z )/a° (10) 



in the perturbation equations. The parameters f2 gas and v are not unrelated. From the statistical 
distribution of the eigenvalues of the deformation tensor it follows that only in 8% of the cases 
there is expansion along all three principal axes®, hence only about 8% of the initial Chaplygin 
gas should fail to condense. Therefore, we expect the initial fraction of uncondensed gas to be 



ttgas Vl-V 2 (ftcDM + ^gas Vl - V 2 )' 1 ~ 0.08. (11) 



gas- 



This equation gives an estimate of the parameter v 2 in terms of Q 

In Fig. ^ we compare the CMB angular power spectrum obtained by implementing Q and 
(JTOJ in a modification of the CMBfast^ code with the WMAP data^. Albeit the result is 
preliminary, and by no means a best fit, it is apparent that the CMB data can be described by 
an evolved mixture of Chaplygin gas and condensate with parameters satisfying Eq. (jllj) . 

In Fig. ^1 we exhibit the baryon power spectrum calculated in the two models: the mixed 
Chaplygin gas and the ACDM with the optimal CMB parameters of Fig. ^p. The model spectra 
have been con volve d with the 2dFGRS window function and their amplitude fitted to the power 
spectrum data^l Again, the data is well described in the range of k in which the window 
function is available. Hence, it may safely be said that unification stands up as a viable scenario. 
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